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BOUNDARY CONCENTRATIONS ON SEGMENTS 


WEIWEI AO, HARDY CHAN, AND JUNCHENG WEI 


Abstract. We consider the following singularly perturbed Neumann problem 

du 

e^Au — u + = 0 , u > 0 in fl, — = 0 on dfi, 

oiy 

where p > 2 and D is a smooth and bounded domain in We construct a new class 
of solutions which consist of large number of spikes concentrating on a segment of the 
boundary which contains a local minimum point of the mean curvature function and has the 
same mean curvature at the end points. We find a continuum limit of ODE systems governing 
the interactions of spikes and show that the mean curvature function acts as friction force. 


1. Introduction and statement of main results 

1.1. Introduction and Main Results. In this paper, we establish new concentration 
phenomena for the following singnlarly pertnrbed elliptic problem 

{ —e^Au + u — = 0 in fl, 

u > 0 in fl, (1.1) 

= 0 on dfl, 

where O is a smooth bonnded domain in with its nnit onter normal u, and the exponent 
p is greater than 2, and e > 0 is a small parameter. We prove the existence of solntions 
concentrating on segment of dfl. 

This eqnation is known as the stationary eqnation of the nonlinear Schrodinger eqnation: 

~ dll) h ^ 

ih—=-—Aijj+ Vijj - ^\ijj\P~‘^ijj (1.2) 

where h is the Plank constant, V is the potential and 7, m are positive constants. Then 
standing waves of (1.2) can be fonnd setting ip = e*®*/^n(T) where E is a constant and the 
real fnnction v satishes the elliptic eqnation: 

— Ji^Av + Vv = \v\^~‘^v (1.3) 

for some modihed potential V. If we consider h —)■ 0, the above equation becomes a singularly 
perturbed one. It can also be viewed as a stationary equation of Keller-Segel system in 
chemotaxis and the Gierer-Meinhardt biological pattern formation system [22], [16]. 

Even though simple-looking, problem (1.1) has a rich and interesting structure of solu¬ 
tions. For the last hfteen years, it has received considerable attention. In particular, various 
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concentration phenomena exhibited by the solutions of (1.1) seem both mathematically in¬ 
triguing and scientihcally interesting. We refer to three survey articles [32], [33] and [38] for 
backgrounds and references. 

In the pioneering papers [34, 35], Ni and Takagi proved the existence of least energy 
solutions to (1.1), that is, a solution with minimal energy. Furthermore, they showed in 
[34, 35] that, for each e > 0 sufficiently small, has a spike at the most curved part of the 
boundary, i.e., the region where the mean curvature attains maximum value. 

Since the publication of [35], problem (1.1) has received a great deal of attention and 
signihcant progress has been made. For spike solutions, solutions with multiple boundary 
spikes as well as multiple interior spikes have been established. (See [2, 3], [5], [6], [10]- 
[13], [17]-[20], [23]-[24], [36], [37], [39]-[40] and the references therein.) In particular, it was 
established in Gui and Wei [18] that for any two given integers k > 0,1 > 0 and k + I > 0, 
problem (1.1) has a solution with exactly k interior spikes and I boundary spikes for every e 
sufficiently small. Furthermore, Lin, Ni and Wei [25] showed that there are at least (^h^^^djv 
number of interior spikes. The hrst author and the third author and Zeng [4] obtained the 
optimal bound of number of interior spikes ^ for general smooth domain in 

A general principle is that for interior spike solutions, the distance function from the 
boundary dkl plays an important role, while for the boundary spike solutions, the mean 
curvature function of the boundary plays an important role. 

It seems natural to ask if problem (1.1) has solutions which “concentrate” on higher dimen¬ 
sional sets, e.g. curves, or surfaces. In this regard, we mention that it has been conjectured 
for a long time that problem (1-1) actually possesses solutions which have m—dimensional 
concentration sets for every 0 < m < N — 1. (See e.g. [32].) Progress in this direction, 
although still limited, has also been made in [1, 14, 15, 26, 28, 29, 30, 31, 41, 42]. For so¬ 
lutions concentration on interior higher dimensional sets, we mention the results in [41, 42] 
where the third author and Yang constructed solutions concentrating on line segment in the 
interior of the domain G. For boundary concentration solutions, we mentioned the results 
of Malchiodi and Montenegro [28]-[30] on the existence of solutions concentrating on the 
whole boundary or arbitrary components of dQ when G C , and solutions concentrating 
on closed geodesics of dQ when G C and also the results of Mahmoudi and Malchiodi [14] 
on the existence of solutions concentrating on the k submanifold of dQ G provided that 
the sequence e satishes some gap condition. The latter condition is called resonance. 

In all the above mentioned papers, for higher dimensional boundary concentration solu¬ 
tions, the concentration sets are either the whole boundary or submanifold of the boundary. 
A natural question is: 

Does problem (1.1) have solutions concentrating on a segment of the boundary? 

In this paper, we give an affirmative answer to the above question. We construct solutions 
concentrating on an open segment 7 of the boundary dkl C provided that 7 satishes the 
following condition: 
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{Hi). Let 7 = 7 ([ 0 ,&]) be a segment of the boundary parametrized by arc length, and 
H{P) be the curvature function of dVL at P. Denote by 

H'i'^is)) = ^i/(7(s)), i/"(7(s)) = ^H{-f{s)). 

Assume that H"{'y{s)) > cq > 0 for all s G [0, b], and H'{j{s))ds = 0. 

Remark 1. From assumption {Hi), one can see that 7 must contain a non-degenerate local 
minimum point of the curvature H. The other condition H'{'y{s))ds = 0 is equivalent to 
H{j{0)) = H{'y{b)), i.e. the curvature at the two end points of 7 must be the same. 

Our main result in this paper can be stated as follows: 

Theorem 1.1. Assume that 7 satisfies {Hi), then there exists Eq > 0 such that for e < Eq, 
there exists boundary spike solution to (1.1) concentrating on 7 . 

Remark 2. In [19], Gui, Wei and Winter proved the existence of multiple spike solutions 
concentrating at the local minimum point of the curvature function H{P). In this paper, we 
have proved the existence of spike solutions concentrating on the segment which contains a 
local minimum of H{P). Theorem 1.1 implies that we can extend their result to a segment 
containing a local minimum point of H. 

1.2. Description of the construction. The solution we construct consists of large number 
( 0 (^j^)) of spikes distributed along the segment 7 whose inter distance is sufficiently small 
((9(eln£)). At the hrst glance one may discard such kind of solutions as there seems to 
be no balancing force at the end points of the segment. In the following we show that the 
derivative of the mean curvature function acts as frictional force. This new phenomena was 
hrst discovered in the construction of CMC surfaces by Butscher and Mazzeo [7]. We will 
comment more on this in later section. 

In this subsection, we will briehy describe the solutions to be constructed later and will 
give the main idea in the procedure of the construction. 

To be more specihc, let w be the unique solution of the following equation: 

{ Aw — w + = 0 in 

w > 0, t(;(0) = max 2 ^gR 2 w{y), (1.4) 

re —)■ 0 as |i/| —)■ 00 . 

It is well-known (see [21]) that w is radial, i.e., w = w{r) and w'{r) < 0 and has the 
following asymptotic behaviour: 

»'(9)=c„j9rAie-l>'l(l + o(l)) (1,5) 

and 

'>^'{y) = -(1 + o(l))w(?/) as \y\ cx). ( 1 . 6 ) 

For P G we set 


fig — Ez G fi}, fi£,p — -P G fi}. 
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and 


where Vn^ pW{z 


Vwp{z) = Vn„w{z - ),wp{z) = w(z -), z e fie 

e e 

— j) is defined to be the unique solution of 
P dll 

Au — u + {w{- - )Y = 0 in fls^p, — = 0 on dfle,p- 


(1.7) 


We will put large number of boundary spikes along 7 . Let the location of spikes be 
( 7 ('Si), • • ■ , 7 (sfc))- We define 






w[z 


i=l 


e 


) 


to be a approximate solution. A natural and central question is how to choose Sj such that 
U is indeed a good approximation. By formal calculation, one has the following energy 
expansion for the energy functional: 


J{U) 


i=\ 


7i . 7 (^ 1 ) -7(gi) ^ 
2 ^ e ’ 


+ o(e) 


where 70 ,7i are positive constants. One needs to find a critical point (si,--- , s^) of J in 
order to get a solution of (1.1), i.e. = 0 for all i. The main point in this paper is to 

exploit the contribution of H'{'y{s)) in The novelty of this paper is the new method of 
constructing balance approximate spike solutions, i.e. the configuration space {(si, • • • , s^)}, 
such that is almost 0 . 

USl 

It turns out that the number of spikes and their positions are determined by some nonlinear 
equations which involves the interaction of spikes and also the effect of the boundary curva¬ 
ture. To explain this, we need to introduce the interaction function \l/(s) which is defined for 
all s G M by 

^(s) = - / w{y - (s, 0))pwP~^—dy. 

JrI oyi 

It turns out that 7 (si) are determined by the following non-linear system: 

' Vl>(h 2 ^)+£ 2 ^'( 7 (s,)) =0(£3)^ 

+ £^H'{^{s 2 )) = 0{eY, 

< ( 1 . 8 ) 

= 0{eY, 

^ = 0{eY, 

and the number of spikes depending on e: is given hj k = ke = + 1 - 

In general, the above nonlinear system is difficult to solve. Our new idea is to consider 
this non-linear system as a discretization of its continuum limiting ODE system (as the step 
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size tends to 0): 


1 

r f = 

1 | = //'( 7 (^(t))), 0 <t< 6 „ 

1 p(0) = 0, p{he) = Pb, 

[ x\K) = -i-^^~\e^H'{j{x{h)))) 

(1.9) 

where is the inverse function of T, and = (/Cg — l)/i = 6 + 0 (/i) and p;, < 0 is a small 

constant depending on e. The above overdetermined ODE is solvable under the assumption 
of the segment 7 in {Hi). 

To describe the configuration space of 7 ( 5 *), we solve the ODE system (1.9) first and denote 
the solution as x{t). Then we define 

c 

H - x{ 2 ) for z - I,--- 1 

(1.10) 

and 

sl = sl_i + e^-\e^H\i^^Pb)) 

(1.11) 

where 

fj = (z — 1) eIne , z > 1 

(1.12) 

Letting yi G M, we define 

Si = s° + yi, for z = 1, • • • ,k, 

(1.13) 


and yi satisfies 


\yi\ < C'|£:ln(-ln£)|, 

K'Si+l ^i) 'Si—1)| ^ 


-- 




^+l 


e 


)} 


(1.14) 


for i = 2, • • • ,k — 1 for some constant C > 0 large. 


With these notations, we can define the configuration space of (.si, • • • , Sk) by 

Afc = {(si,-- - , Sfc) G M^|sj is defined by (1.13) and satisfies (1.14)} (1-15) 


Moreover, from the analysis of the ODE (1.9) in Section 6, one can get that 

9 • — ■ 1 CP 

|si - Si-i\ > (1 + o(l))|£lne|, w{— - —) < - —r (1-16) 

e |ln£| 

for i = 2, • • • , fc and 

|si - Si_i| = 2(1 + o(l))|£lne|| 

for = 1, A:. 

We will prove Theorem 1.1 by showing the following result: 


(1.17) 
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Theorem 1.2. Let j be a segment of dQ satisfying the assumption {Hi). Then there ex¬ 
ists £0 such that for e < Eq, there exists positive number k = = [| 7 ]|-^] o-nd k points 

( 7 (si), • • • , 7(sfc)) on 7 , where (si, • • • , G such that there exists a solution to Prob¬ 
lem (1.1) and Ue has the following form: 


k 


uAx) = y V, 


i=l 


7(si) 


w 


.x - j(Si), 


+ 0 ( 1 ) 


where o(l) —)■ 0 as £ ^ 0 uniformly. 


(1.18) 


Remark 3. The motivation of our construction comes from the study of the constant mean 
curvature surface. In [7], Butscher and Mazzeo constructed CMC surface condensing to a 
geodesic segments by connecting large number (0(^)) of spheres of radius r distributing 
along the geodesic segment. In their paper, they require the symmetry condition on the 
geodesic segment. In our main theorem 1.1, if we further require that hi is symmetric, it is 
easy to see that (Hi) can always be satished near the non-degenerate minimum point of the 
curvature H{'y{s)). Since Theorem 1.1 can deal with more general segment, we believe that 
our idea can be used to construct CMC surface condensing to geodesic segments without the 
symmetry condition. We will discuss this in a forthcoming paper. (A. Butscher announced 
this result in the preprint [9] but the full details have never appeared.) 

1.3. Sketch of the proof of Theorem 1.2. We introduce some notations hrst. Since after 
scaling x = ez, the original problem is 


Au — a -I- = 0 in fie 

M > 0 in fl 

= 0 on an. 


Fixing s = (si, • • • , Sfc) G A*,, we denote by 


P = (Pi,---,P,) = (^,---,^) 


5 5 


and dehne the sum of k spikes as 


U = '^Vn^p.w{z- Pi). 


2=1 

Dehne the operator 

S{u) = Au — u-\-u^. 

We also dehne the following functions as the approximate kernels 

_ - C) 

Zj^ , t 1, , /i. 

UTp, 


(1.19) 

( 1 . 20 ) 

( 1 . 21 ) 

( 1 . 22 ) 

(1.23) 


Using U as the approximate solution, and performing the Lyapunov-Schmidt reduction, we 
can show that there exists £o such that for £ < £o, we can hnd a of the following projected 
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problem: 

k „ 

S{U+ 'ip) = '^CiZi, / = 0, z = 1, • • • , A:, (1-24) 

i=l 

where Cj are constants depending on the form of 'ip, Zi. 

Next, we need to solve the reduced problem 

Ci = 0, i = ,k (1.25) 

by adjusting the points in A^. 

There are two main difficulties in solving the reduced problem. First we need to control 
the error projection caused by ijj. In order to control this projection, we need to work in 
a weighted norm, which estimates -0 locally (see Section 4), and also we need a further 
decomposition of 0 which is given in Section 4 from where one can see why we dehne the 
conhguration space of Si in (1.15). Second, we need to solve a non-linear system of the form 
(1.8), for which we use the discretezition of the ODE equation (1.9). 

Finally, the paper is organized as follows. Some preliminary facts and useful estimates are 
explained in Section 2. Section 3 contains the standard Liapunov-Schmidt reduction process: 
we study the linearized projected problem in refsec4.1 hrst and then we solve a non-linear 
projected problem in 3.2. In Section 4 we obtain a further asymptotic behavior of 0 in e. 
In Section 5, we derive the reduced nonlinear system of algebraic equations for the location. 
Section 6 is devoted to solve the nonlinear system. 

Acknowledgments. The research of J. Wei is partially supported by NSERC of Canada. 

2. Technical Analysis 


In this section, we introduce a projection and derive some useful estimates. Throughout 
this paper, we shall use the letter C to denote a generic positive constant which may vary 
from term to term. By the following rescaling 



X = ez, z E fie ■= G D}, 

(2.1) 

and equation (1.1) becomes 

f Au — u + = 0 , in De 

\ Ip = 0 on dfle. 

(2.2) 

We denote by = {(z/i. 

y2)\y2 > 0}. Let w be the unique solution of 



( Aw — tc -|- = 0 in 



< re > 0, r(;(0) = maXj^g^2 w{y), 

[ w{y) — )■ 0 as \y\ —oo. 

(2.3) 


Let p G dQ. We can dehne a diffeomorphism straightening the boundary. We may assume 
that the inward normal to dQ at p is pointing in the direction of the positive X 2 axis. Denote 
B'{R) = {|xi| < R}, and Di = QnB{p,R) = {(xi,X2) G B{P,R)\x2 - P2 > p{xi -pi)} 
where B{p, R) = {x E R‘^\\x — p\ < R}- Then since dQ is smooth, we can End a constant R 
such that dfl can be represented by the graph of a smooth function pp : B'{R) —)■ M where 
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Pp(0) = 0, and /9p(0) = 0. From now on, we omit the use of p in pp and write p instead if 
this can be done without confusion. So near p, dVt can be represented by (xi,p(xi)). The 
curvature of dVL at p is H{p) = p''{0). After scaling, we know that near P = |, can be 
represented by {zi,e~^p{ezi)), where {zi,Z 2 ) = e~^{xi,X 2 )- By Taylor’s expansion, we have 
the following: 


e V(ezi) = ip"(0)£2:f + ip<‘‘>(0)eVJ + Ois'-’zf). 


(2.4) 


Recall that for a smooth bounded domain U, the projection Pjj of H‘^{U) onto {u G 
H'^{U)\^ = 0 at dU} is dehned as follows: For v G H‘^{U), let Puv be the unique solution 
of the boundary value problem: 


Au — u + yP = 0, in U, 
g = 0 on dU. 

Let hp{z) = w{z — P) — Pn^pw{z — P), then hp satishes 

Ahp{z) — hp{z) = 0, in fig, 

dhp 
du 


^ on aflg. 


= 


For z G fli^g, for P = (Pi, P 2 ), set now 


2/1 = A 
2/2 = 2^2 


Pi, 

P2-£-V(£(a-Pi)). 


(2.5) 


( 2 . 6 ) 


(2.7) 


Under this transformation, the Laplace operator and the boundary derivative operator 
become 

A, = Ay + p{ezifdy^y^ - 2p\ezi)dy^y^ - ep"{ezi)dy^, 
d , 2 , 


(1 + p'{eziY) 2 — = p'{ezi)dy^ - (1 + p' {^zi))d, 


y2- 


Let uB) be the unique solution of 


An — n = 0, in 


(2.8) 


where w' is the radial derivative of w, i.e. w' = Wr{r), and r = \z — P\. 
Let be the unique solution of 


An — n — 2p''{0)yi 


= 0 in 




W 2 = 

Let n^^^ be the unique solution of 

An — n = 0, in 




(2.9) 


fe = hF"‘'(0)p?. on aiR^. 

Note that are even functions in yi and n*-^^ is odd function in yi. Moreover, it is 

easy to see that \vi{y)\ < for any 0 < p < 1. Let x(^) be a smooth cut-off function. 
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such that xi^) = 1, a: G B{0,RQe\ lne|), and x(a;) = 0 for x G i?(0,2i?oe| lne|)‘^ for Rq large 
enough, and Xei^) = xi^^) for ^ ^ fo this case, one has ^(-Rol fo^l) = Set 

hp{z) = -{evi{y) + e^{v 2 {y) + v:i{y)))xe{z - P) + e^^p{z), z G (2.10) 

Then we have the following estimate: 


Proposition 2.1. 

||^(^)||Hl(Oe) < C"- (2-11) 

Proposition 2.1 was proved in [39] by Taylor expansion and a rigorous estimate for the 
reminder using estimates for elliptic partial differential equations. Moreover, one can checked 
that |^(^)| < for some 0 < /i < 1. 

Similarly we know from [39] that 


Proposition 2.2. 


dw dVnW 


](z - P) = er]{y)xe{z - P) + £^Vi{z), z G 


drp drp 

where rj is the unique solution of the following equation: 

At] — r] = 0 in R^, 

w, = -Kfp - - ii/CDw SRl 


Moreover, 


Ihillr/pOe) ^ 


( 2 . 12 ) 


(2.13) 

(2.14) 


One can observe that r]{y) is an odd function in yi. It can be seen that \r]i{y)\ < 
for some 0 < /i < 1. 

Finally, let 

Lo = A — 1 + pw^~^{z). (2.15) 

We have 


Lemma 2.1. 

(iin 

Ker{Lo) 0 iP^(R+) = span{ — }, 

oyi 

where = {n G ^ = 0 on 

Proof. See Lemma 4.2 in [34]. 

Next we state a useful lemma we will always use: 


Lemma 2.2. If |gi — q 2 \ « | 5 i|, we have the following estimate: 
f div 

/ PU!{yY~\w{y - qiCi) + w{y + q 2 ei))^dy = 0{\qi - q 2 \w{\qi\)) 
JRI dyi 


(2.16) 


□ 


(2.17) 


as |gi| —)■ cxo. 
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Proof. By the oddness of in ?/i, one has 


f dvu 

/ Pw{yY~\u!{y - qiei) + w{y + q2ei))—dy 
jRl dyi 

= [ Pw{yY~^iwiy - qiei) - w{y - q2ei))^dy 
JRI oyi 

f dxv 

= / Pw{yY~^\-^\ 0 {w\y - qiei)\qi-q2\)dy 
JR2 oyi 

= 0 {\qi - q2\)w{\qi\). 


□ 


Remark 4. In the following sections, we will denote by y^ = {y\, yf) to be the transformation 
defined by (2.7) centered at the point Pi and be the corresponding solutions in the 
expansion of hp^. 


3. Liapunov-Schmidt Reduction 


In this section, we reduce problem (2.2) to finite dimension by the Liapunov-Schmidt 
reduction method. We first introduce some notations. Let P[ff{Qs) be the Hilbert space 
defined by 

f)q I 

= {ue HYQYItt = 0 on (3.1) 

ou 

Define 

S{u) = Au — u + u‘^ (3.2) 

for u G P[jf{Qe)- The solving equation (2.2) is equivalent to 

S{u) =0,ue (3-3) 


To this end, we hrst study the linearized operator 

k 

LsiY) ■= Af; - fj + pC^Vn,^p,w{z - Pi)y-^'ilj, 


2 = 1 


and dehne the approximate kernels to be 

dVn„w{z - Pi) 


Z, = 


dp 


Pi 


for i = 1, • • • ,k. 


(3.4) 


3.1. Linear projected problem. We first develop a solvability theory for the linear pro¬ 
jected problem: 

( = h + Y!l=iCiZi, 

{ L fjZidz = 0,z = 1, • • • 

1 e Hiyiy 


(3.5) 
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Given 0 < /x < 1, consider the norms 

= (3.6) 

z&ile j 

where Pi G with dehned in (1.15). 

The proof of the following Proposition on linearized operator, which we postpone to the 
appendix, is standard. 

Proposition 3.1. There exist positive numbers p G (0,1), Sq and C, such that for all e < Eq, 
and for any points Pj, j = 1,... ,k given by (1.15), there is a unique solution cf) to problem 
(3.5). Furthermore 

||^IU<q|h||.. (3.7) 

In the following, if is the uniqne solution given by Proposition 3.1, we set 

f; = A{h). (3.8) 

Estimate (3.7) implies 

\\A{h)\U < C\\h\U. (3.9) 


3.2. The Non-linear Projected Problem. We are now in the position to solve the equa¬ 
tion: 

r L,^^; = E + N{^|;) + Zl,c^Zi, 

< ifZi = 0 for i = 1, • • • , fc, (3.10) 

{ G 

where E is the error of the approximate solution U: 

k k 

E = ( 3 . 11 ) 

2=1 2=1 

k 

+( 5 ;Po,.,.w(.'-P.))’’, 

2 = 1 

and Nlfj) is the nonlinear term: 

k k 

NW = - m + i>r - (Y,rn.,,w(z - Pi))” ( 3 . 12 ) 

2=1 2=1 

k 

2=1 

We have the validity of the following result: 

Proposition 3 . 2 . There exist p G (0,1), and positive numbers Eq, C, such that for all e < Eq, 
for any points Pj, j = 1,... ,k given by ( 1.15), there is a unique solution {fj, Cj) to problem 
(3.10). This solution depends continuously on the parameters of the construction (namely 
Pj, j = 1,... ,k) and furthermore 


IN’II. < Ce. 


( 3 . 13 ) 
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Proof. The proof relies on the contraction mapping theorem in the || ■ ||*-norm above intro¬ 
duced. Observe that solves (3.10) if and only if 

= A{E + (3.14) 

where A is the operator introduced in (3.8). In other words, solves (3.10) if and only if 
is a hxed point for the operator 

Ti'iP) := A{E + N{f:)) . 


Given G > 0 large enough, dehne 


B = //^(a) : < Ce, [ fjZi = 0} 

J Qg 


We will prove that T is a contraction mapping from B in itself. 

To do so, we claim that 

\\E\\,<Ce (3.15) 

and 

\\Nmu<c\mi ( 3 . 16 ) 

for some hxed function C independent of e, as e ^ 0. We postpone the proof of the estimates 
above to the end of the proof of this Proposition. Assuming the validity of (3.15) and (3.16) 
and taking into account (3.9), we have for any E B 


l|r(V>)ll.< c[||i5 + iV(^,)||.]<c(£ + £^) 

< Ci£ 


for a proper choice of Ci in the dehnition of B. Take now fji and '02 in B. Then it is 
straightforward to show that 

||r(0i)-r(02)ll* <q|iv(0i)-iv(02)||* 

< C [||0l||* + ||02||*] 1101 - 02||* 

< O ( l )||01 - 02 ||* 

This means that T is a contraction mapping from B into itself. 

To conclude the proof of this Proposition we are left to show the validity of (3.15) and 
(3.16). We start with (3.15). 

Fix j E {1,k} and consider the region \z — Pj\ < ™™hG-A-d’IA--Pj+i|} ^ region 

the error E, whose dehnition is in (3.12), can be estimated in the following way 

IBWI < CwO-^z - Pi)lY, + 

Pi^Pj i 

< C{e + 

if we choose p small enough such that p — n> 2. 


(3.17) 
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Uonsider now the region \z — Fj\ > -— 2 — 

E, we get in the region under consideration 


, for all j. From the definition of 


\E(z)\ < C 


E + (E - ■P*))" - E - P<)' 


i=l 




< 


c^E' 


-lJ,\z-Pi\ 


(3,18) 


From (3.17) and (3.18) we get (3.15). 

We now prove (3.16). Let E B. Then 


k k 

iAf(v>)i < \{(j2Pa..r«u-m+i>r-(E,p^Mz-Pi)r 

i=l i=l 

k 

-p(J2p«:pM- - p.)rvi < cii\ 

i=l 

Thus we have 
This gives (3.16). 

A direct consequence of the fixed point characterization of ijj given above together with 
the fact that the error term E depends continuously (in the *-norm) on the parameters Pj, 
j = 1,... ,k is that the map (Pi, • • • , Pk) into the space Cipl;.) is continuous (in the 

*-norm). This concludes the proof of the Proposition. 

□ 


4. Further expansion of the error 


For later purpose the asymptotic behaviour of the function -0 as e ^ 0. This is needed to 
compute the neighboring interactions. 

Before we state the result, we first consider the following equation: 


— 0 + pwiyY ^(j) = h + 
^ = 0 on dRl, 

^^dy = 0 

oyi ^ 


in 






(4.1) 


where d = 


J^2 

r Swl-). ■ We consider the above equation in the space where 

APat/i) 


< + 00 , where 


. = supj^gR^ for some 0 < < 1. It is quite standard to show the solvability of 

the above equation and (j) will satisfy the following estimate: 


Ill’ll.. < C||ft|| 


( 4 . 2 ) 








14 


W. AO, H. CHAN, AND J. WEI 


Now we decompose ip as follows: 

Proposition 4.1. 

k 

'P = - Pi)(pi + e'^ipi, 

where 


2=1 




and (pi = pi{y^) is the unique solution of 


Api - pi+ pw{ypP Vi = Hi + in 

= 0 on dM.1, 
oy2 +’ 


El. 


(4.3) 

(4.4) 

(4.5) 


ful = 0 

where di is defined such that the right hand side of the above equation is orthogonal to ^ 
in Lfi norm, and 


Hi = -pw{y'-Y Yw{y" - —^ei) + w{Y - —^ei) + (4.6) 


for i = 2, - ■ ■ ,k — 1 and 


Hi = -pw{yY^ ^[w{y^ - — —^ei) + 


and 


and we denote 


Hk = -pw{y^r-Yw{y^ - ^ 


-ei) + eV ]> 


are the solutions obtained in Section 1.2 centered at the point Pi. 
Proof. First by the definition of di, 

d.= f H,^dy 

Jr^. dyi 


(4.7) 

(4.8) 

(4.9) 


(4.10) 


Then from Lemma 2.2, and the evenness of with respect to y\, and the dehnition of the 
conhguration space (1.15), we know that for i = 2, - ■ ■ ,k — 1 


\di\ < Ce V-Si+i - Si\ - |si - Si_i| min{w;(^^^——fhJ. 


)} < CV, (4.11) 


and for i = l,k, 

|di| = 0{w{— -^)) = 0{eY and \dk\ = 0{w{- 

Moreover, from (4.1), we have the following estimate: 

\\Pi II *4= < Ce if p > 2 + ^1. 


-)) = 0(e). (4,12) 


(4.13) 
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Our strategy to estimate "01 is to decompose 0i into three parts and show that each of 
them is bounded in || • ||* as e ^ 0. We write 0i as 


where 0ii satishes 


Dehne 0i2 by 


01 = 011 + 012 + 013, 


A 011 - 0 n = 0 , in Oe 


(4.14) 


(4,15) 


'^’12 = 3 


2 = 1 


and Si is determined by 


M{si) = 


« k 

- / iY] Xs{z - -Pi)0i + 


Finally dehne 0i3 to be the solution of the following equation: 

^0013) = ^^00 - EJ =1 Xe{z - - £0011 + 012)) m 

^ = 0 on ao, 

/ 013^4dz = 0. 


(4.16) 


(4.17) 


(4.18) 


Next we will estimate 0ii, 0i2, 0i3 term by term. First we estimate gi^ = 


By direct calculation 


du 


9le 


1 1 ^ ^Xe{z - Pi)-, 

^ Y^Xeiz -P^)^ + 0-EiT-) 


2 = 1 
k 


ee 


0 ^Xe{z - Pi) 

-2 

2 = 1 


dix 

k 


dix 


+ 0(0) 


(9('£-2g-(/4l-/4)^o|lne|^ g-/4|2-Pi| 


2=1 


< 






2 = 1 


if we choose ni > g and the cutoff function in such a way that (pi — /x)i?o > 1- In the 
above estimate, we use the dehnition of 0* and the Neumann boundary satished by it and 
the dehnition of the cut-oh function y. Thus we have that ||5'ie|l* < C, therefore, there exists 
constant C > 0, such that 

||0ii||4<C'. (4.19) 
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By the definition of ' 0 i 2 , (f)i and the estimate on one can obtain that 


/ Cyi Xe{z - Pj)(t)j + 

= / Xe{z - Pi)4>iZidz + ^ Xe{z - Pj)4>jZidz 

+ 0{e^) 

In order to estimate the above term, we first consider a general fnnction which is the 
solution of the following equation: 

A(f) — (f) + pw{yy~^(j) 

= -pw{yy-^{w{y - Qiei) + w{y + gsCi) + + d^^ 

' ^ = 0 on dM.1, 

. fui = 0 

We can decompose it as 


dyi ^+’ 


(4.20) 


where 


and 


= 


" A(f)^ — (j)^ + pw{yy 

= —pw{yy~^{w{y — qiei) + w{y + qiei) + 

g^ = 0on9M^, 

. y = 0 


" AcJ)"^ — (j)"^ + pw{yy ^(jZ 

= -pw{yy-^{w{y + 9261 ) - w{y + gid)) + d 2 ^^ 


dy, 


(4.21) 


(4.22) 


It = 0 

fmi = 0 


a— m Kl, 
ayi +’ 


(4.23) 


.^j]5.q_ ' dyi 

where di are defined such that the right hand sides of the above equations are orthogonal to 
^ in norm It is easy to see that (j)^ is even in yi and by Lemma 2.2, cfZ satisfies 

< Cw{qi)\qi - ^ 2 !, 

if \qi - q 2 \ « |gi| and |gi| -)■ cx). 

Using the above estimates, we can decompose (pi as 


— Pi.l + 


and pi^i is even in y\ and 
||0i,2||** < C\\ 


‘i,l “T (Pi,2 


^ i —1 I I ^ 2+1 11 f ^ 1 \ „ '^ 2+1 ^ 


min{t(j(- 




(4.24) 
)} < Ce\ (4.25) 


Then by the above estimate and the decomposition in Proposition 2.2, we have 
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Pi)(j)iZidz = 0{e^), 


(4.26) 


and similar to the decomposition of (j)i, one can also decompose + as an even function 
of yl and an O(e^) function, so we get that 


/ Xe{z - Pj)(j)jZidz 


= f {(pi-i + dy + 0{e^) 

jRl dyi 

= 0(8^) 


Moreover, since |si — S 2 I = 2(1 + o(l))|e lne:| and |sfc_i — Sfc| = 2(1 + o(l)) |e lne|, one can get 
that 



P2)(t)2Zidz = 0 (£^), [ Xeiz 
J fie 


Pk-i)(t)k-iZkdz = O(e^). 


(4.27) 


Thus we have 


Si| < Ce^. 


(4.28) 


Next we estimate '0i3- Denote by 


k 

fe = L.ii; -^Xeiz - Pi)(pi - 8^ + iJu)). 
i=l 


Claim: 


II/. 


< Ce/ 


(4.29) 


Proof of the Claim: 
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By the definition of f^, we have 

k 

f^{z) = Lei'll; - ^ Xe{z - Pi)4>i - £^('011 - ' 012 )) 


i=l 


= E + iV(0) + ^ L,{xe{z - PM - + M 

i i 

= (E - -p*))' - E®*' - P‘)' + + E<=‘^- 

i i i 

- X] Xe{z - Pi){Mi - PneW{z - Pi)Y~^ + O(e))0i) 

i i 

+ Y^{2VM {Xeiz - Pi)) + Mxeiz - Pi)) - e^Leilpu + '012) 

1 

= (E - p*))” - E - py ++ E '=‘^- 

i i i 

- Y, Xei~ - ^’^)(p((E P^..pM^ - P*))”"' - ^iv - P)”"')* 

2 i 

-pw{y - Pi)P~^{w{y - Pi_i) + w{y - F^+i) + evuiy)) + diZi) 

+ ^O(£)0i + M{Xs{z - Pi)) + (f)iAXe{z - Pi)) - e'^Lei'lpll +012). 

i i 

From the definition and estimates of 0*, 0ii, 0i2, x, and the configuration space, we know 
that |cj| = 0(£^), so 

11/00 < C '£0 

By the a priori estimate, we know that 

II013II* < C, 

thus we have 

ll^ilUEa 

We thus finish the proof. 

□ 


Given points Pj defined by (1.15), Proposition 3.2 guarantees the existence (and gives 
estimates) of a unique solution 0, c*, z = 1,..., /c, to Problem (3.10). It is clear then that the 
function n = f/ + 0 is an exact solution to our problem (1.1), with the required properties 
stated in Theorem 1.2 if we show that there exists a configuration for the points Pj that 
gives all the constants Cj in (3.10) equal to zero. In order to do so we first need to find the 
correct conditions on the points to get Cj = 0. This condition is naturally given by projecting 
in LP‘{Vti;) the equation in (3.10) into the space spanned by Zj, namely by multiplying the 
equation in (3.10) by Zi and integrate all over Gg. We will do it in details in the next final 
Section. 
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5. The Reduced Problem 

In this section, we keep the notations and assnmptions in the previous sections. As ex¬ 
plained in the previous section, we have obtained a solution u = p ~ + 

XliLi Xe{z — Pi)(f)i -|- e^ipi of the following equation 


r d^u-u + vP = Ci-Z’i in hie 
\ = 0 on (?hle 


In this section, we are going to solve c* = 0 for all i by adjusting the position of the spikes, 
i.e. Pi- First, multiplying the above equation (5.1) by Zj, i = 1, • • • , A: and integrating over 
Re, we get that 


■ Cl ■ 



C2 

= 

/^^(AM-n + nP)Z2 

_ _ 




(5.2) 


Recall that M is invertible, so Cj = 0, i = 1, • • • , /c is reduced to solve the following system: 


- u + uP)Zi 
-u + uP)Z 2 

L_(A«-m + u'’)Zs 


(5.3) 


We have the following estimates: 

Lemma 5.1. Under the assumption of Proposition 3.2, for e small enough, the following 
expansion holds: 



u -|- u^)Zidz 


and for i = 2, - ■ ■ ,k — 1, 


- e^iz 2 H'{j{si)) + 0{e^), 
e 


(5.4) 



u + u’')Z,dz = —^tl) - —5iti) - eT 2 ff'( 7 (sj)) + O(e^) 

£ £ 


(5.5) 


and 



(An — u + u^)Zkdz = 


- £^i'2H'{'y{sk)) + O(e^). 


where 1^2 > 0 is a constant defined in (5.12). 


(5.6) 
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Proof. First, by direct calculation, one can get the following expansion: 

Am — u + 

k k 

A{U + Xe{z - -{U + Y, Xe{z - Pi)4>i) + {U + Y Xe{z - 

i=l i=l 

k 

+ £^(A?/;i - + p{U + Y Xe{z - Pi)(piY~^'i/ji) 

i=l 

k k 

(u+Yl -(^+Y1 xe{z - md' 


i=l 
k 

ViU + Y Xe{.Z - Pi) 0 i)^"VVl 

i=l 

:= Ii + I2 + h- 


2=1 


Next we calculate Ii to Is term by term. First from the estimate on fji in (4.4) 

f l2Zi = e^[ - -01 + p{U + Y\ Xejz - Pi)(t)iY~^'il)i)Zi 

J Qg- J fig ^ 

-pw{z - - p^m^-^zpf. 


= £ 


= £^ / p{p - 1 )m;(2 ; - Pi} 


dr 

,_2dw{z - Pi) 


dr 




j¥=i 


dw{z — Pj) 


+ 0{e))dz 


Moreover 


= O(e^). 


hZi = 


\u + Y. - ('^ + E - a)*)' 


(5.7) 


(6.8) 


-p(U + Xe(z - Pi)4'i)'’ “eVl 

i 

< c[ £Vin^*i = o(£='). 

J fig 

Next from the equation satished by (pi and the dehnition of the cutoff function y, we get that 

f hZ, = f (AU-U + U^Zi + Y f Xe{z-Pj){A(Pj-^j+pUP-^fj)Z, 

J fig J fig j J fig 

+ [(tf + XeO - pm- -U-- pU--' Y.*.(z- B)*lZi + 0 (e=>) 

i i 

= [ (AU-U + UP)Zi + h^ + h2 + 0{e^). (5.9) 
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Similar to the estimate in (4.26), using the equation satished by (pi, we have for i = 2, - ■ ■ , k—1 

+ pU^~^(pj)Zi 

= [ Xe{z - Pi){Ap)j - (pj + pUP~^(pj)Zi 
J £7e 

+ + pU^~^(pj)Zi 




; J fie 


= / p([/^-^-w^')P>^Zi + 0(e^) 

J Qg 

+ / Xiz - Pj)iA(pj - (pj+pUP-^(pj)Zi + 0{e^) 




= / P{p-^)wf {Wi+i + Wi-i)(piZi 

J £7e 

+ E / Xe{z-Pj){A(P,-(P,+pUP-^(Pj)Z, + 0{e^) 




= 0{e\ 


I ^i—1 I I ^i+1 I 


min{i(;( 


^i—l\ '^ 2+1 

- , w; - 


)}) + 0 (£=^ 


= 0 ( 6 =^) 

and similarly we can always decompose 


Xe{z - Pj)(pj = eiJij + 0{e^) 

j=i 

where ipij is a function even in y\, and by Proposition 2.2, we have 

dw{y^) 


Z^ = 


where pi is odd in y\. Thus we have 


dyi 


+ ^Pi + 0{e ) 


Ii 2 <C j p{p - 2)wf ‘^(Y Xe{z - Pj)(pjfZidz + 0{e^) 

Jfle j = i 

< Ce\ 

For the case i = l,k, recall that ) = O(e^), one can also get that 


Ji1+Ji2 = 0(£3). 


Thus we have the following: 


hZ,dz= / {AU-U + U’’)Zi + 0(s^). 


(5.10) 
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Next for i = 2, • • • ,k — 1 


/ {AU-U + UP)Zi 

'Os 

i 

{w{z - Pi) + ^ + "y®) 





>n, 


+E - Pi) +o(£'’))’' - E - ■p*)’' 






= f pw{z - Pi)^ ^ ^ ) 

iOe ^ 

+w{z - Pi_i) + w;(2r - Pi+i)^ 

= f pwiy)iw{y-— — -ei) + w{y -—— 


dyi 


-2 / ^....„ap-1^(2/).,(3) , 


+£ / pw{yy 


dyi 


-vr + 0{e^) 


Similarly, one has for i = l,k, 


{AU -U + UP)Zi 


+r 


S 2 -S 1 .dw{y) 

pw[y)w[y -ei)—- 

£ dyi 

R 2 


and 


/ (A[/-17 + 17?>)Zfc 

Jfle 

f , I , Sk-i-Sk .dw{y) 

= l,My)'ii{y - —^i)^ 

+£' f pw{yY^'-‘^PZ..ifp + 0 (£“). 

JKi oyi 


(5,11) 
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(3) 

Next by the definition of , we can get that 

[ dy = f 

Jm oyi Jm2 




dyi 


dw{y) dv. 


(3) 


( 3 ) d dw{y) 

— V- - -- dy 




JaM.1 dyi dy2 * dy2 dyi 

^ -II 

= -V2P^'^\Pi) = -V 2 H'Yi{Si)) , 

where U 2 = | /]R(j^)^2/i > 0 is a positive constant. 

Recall that the interaction function is defined by 


^(s) = - pw{y - {s,0))w{yY- 


dyi 


Combining (5.9), (5.10),(5.11), (5.12) and (5.13), we know that 

^ iSi —S 2 |^ _2,, I 


hZidz = -^( 1 - 


.\)-e^i.2H'{j{sY) + 0{eY 


and for i = 2, • • • ,k — 1 


[ hz,dz = d/(|^^^|) - vl/(|£^^|) - e‘^U 2 H'{j{sO) + 0{eY 

Jne ^ ^ 

and 

[ hZkdz = - e^U2H'{j{s,)) + 0{eY. 

Jne ^ 

The results follows from (5.7), (5.8) and (5.14)-(5.16). 


From Lemma 5.1, the problem (5.3) is reduced to the following system: 
vl;,(|^|)+£ 2 //'(^(si)) = 0(£3), 

^i(l“l) - ^( 1 ^ 1 ) + e^H\^{s2)) = 0(£3), 

< ; 

^ i (|£^ z £^|) _ ^(1 + e^H'{^{sk-i)) = 0{eY, 

^ _Vl/i(|£^|)+e2i7'(y(s,)) = 0(£3). 


where we denote by 




(5.12) 

(5.13) 

(5.14) 

(5.15) 

(5.16) 

□ 


(5.17) 


By summing up the first i equations, one has 

r ^i( ^i+i-^i ) ^ e^H'{'y{sj)) = 0{ie^) for i = 1, • • • , A: - 1, 

I Y.Ue^H\^{s,)) = 0{keY. 


(5.18) 
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6 . Solving the nonlinear system 


Our aim in the rest of this paper is to find a solution {s*} to the non-linear system (5.18) 
in (1.15). 

Observe that the linearized matrix of the above system at main order is degenerate, thus 
the terms containing H'{j{s)) will play an important role. We will explain how we solve 
system (5.18). The novelty of this paper is to consider the above system as a discretization 
of an ODE. In order to explain this idea, we first introduce some notations. 


Let 

s = G{b) 

be the solution of Ti(s) = b. Since Ti(s) = C'„s“5e“^(l -|- o(l)) as s oo, using this 
asymptotic behaviour of Ti, one has the following: 


G(0 =-(1+ 0(i^^^-4^))ln5, as 6^ 0. (6.1) 

In b 

Then the above reduced system (5.18) is equivalent to the following system: 

f Si+i -Si = eG{- + 0{eH)), for f = 1, • • • , A: - 1 . , 

\ Sk- Sk-i = eG{e^H'{'y{sk)) + 0 {e^k)). 


Let h = —elne be the boot size, if we denote by s* = x{ti) where ti = [i — l)h, then from 
the above system (6.2), 


E’.i + 0(eH)), 

+ 0(A)). 


(6,3) 


In order the solve the above system, we consider the limiting case of the above system, 
i.e. view {'^{x{tj)))h as H'{^{x{t)))dt, and introduce the 

following ODE: 

f = -i^G(j^dW), 

f = H'{^{xm. 

p(0) = 0, p{be) = pb, 

x'{b,) = -^G{e^H'{j{x{b,)))), 
where b, = {k - l)h = [|]h = b + 0 {h). 


(6.4) 


One can see that the above second order ODE has three initial conditions. Besides the two 
end point initial values, there is an extra condition, i.e. the last equation of (6.4), which in 
fact comes from the last equation of (6.2). This ODE with extra initial condition is not always 
solvable. It turns out that this extra condition corresponds to some balancing condition of 
the curvature of the segment 7. In order to solve this ODE, we need assumption {Hi) on 7. 
For this ODE, we have the following existence result: 


Lemma 6.1. Under the assumption (Hi), there exists £0 > 0, such that for every e < Eq, 
there exist pb = Pb{£) < 0, such that the above ODE ( 6 . 4 ) is solvable. Moreover, pb satisfies 






CONCENTRATIONS ON SEGMENTS 


25 


the following asymptotic behaviour: 


ft = -(J/TT.)) + o( ‘"*, 


In e 

Proof. From the asymptotic behaviour of G, we know that the hrst equation of (6.4) is 

dx 
dt 


(6.5) 


■^G(^p(i)) 

In £ In £ 


where 


= (l + 0(hf^))(a,l,i(-p(())+a2) 


1 ln(—hie) 

0-1 = 1 —) 02 = 1 - ] -■ 

In e In e 


Integrating the above equation from be to t, one has 


x{t) - x{be) = 


^ -G{^p{t))dt 


In e In e 


= (l + 0 (M_^))[a 2 (t- 6 ,) + ai f\n{-p{t))dt]. 

Plugging the expression for x{t) into the second equation, 

p'{t) = H\'j{x{be)) + (1 + ^^^^^ ))[a 2 (t - be) + oi f \n{-p{t))dt]). 

In £ Jbe 

By the boundary condition p(0) = 0,p{be) = Pb, we have 

H'{'^{x{be) + (1 + - be) + Oi [ ln{-p{t))dt]))dt = -pb. 


In e 


By Taylor’s expansion, 


H\j{x{be) + (1 + 0 (i^^^j^))[a 2 (f - be) + ai l\n{-p{t))dt])) 

= H\^{x{be) + 02(f -b)+ O(MZM))) 


= H\^{x{be)+a,{t-be))) + 0{ 
So from (6.7) and the above equation, we have 


In e 
ln(—Ine), 
Ine 


H'{j{x{t)))dt = / H'{^{x{be) + a 2 {t - b)))dt + 0{ 


ln(— Ine), 
Ine 


= - a-A)) - + 0(!d-!e£)) 


Ine 


= Pb- 


( 6 . 6 ) 


(6.7) 


( 6 . 8 ) 
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Since by the third boundary condition 



X'(k) = ~G{e^H'{i(x{k)))), 

(6.9) 

one can get that 

Pb = H'{'y{x{bs)))e\ne. 

(6.10) 

We assume that 

pb = {H'{'y{bs)) + pe)e\\ie, 

(6.11) 

then 

x(b)-b + 

(6.12) 

Using (6.12), (6.8) 

is reduced to the following: 



^7 7n^^ ^7 7. , ^'(7(0)) - H'i^ib,)) 

77(7(0)) 77(7(6=))+ P. 

(6.13) 




By the assumption (Hi) 

i/(7(0)) = i/(7(&)), ^ (6.14) 


and 

//"(7 )>Co> 0 , b, = b + 0 {h), 
the above equation is uniquely solvable with 


Pe = 0{ 


ln(— hie) 
Ine 


(6.15) 

(6.16) 


So there exists unique ph = {H'{'^{b^)) + ^”^^ ))£Ine such that (6.4) is solvable, and 


we have 


^,ln(—Ine), ,, , , ^,ln(—Ine), 

^(0) = 0( ). x%) = he + 0{ 


In e 


ln£ 


(6.17) 

□ 


We will use the solution of the ODE to approximate the solution of (6.2). In order to 
obtain a good approximate solution, one need to control the error of 


'•n+i 


j=i 

So we will use the midpoint Riemann sum approximation of integrals which will give us 


U+l 


J 2 H\-f{x{tj)))h- I H'{'y{x{t)))dt = 0 {h‘^). 

j=i 

To be more specihc, we will choose the approximate solution to be the following: 

= x(fi), ti = = 1, • • • , A: - 1, 


(6.18) 


(6.19) 
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and 

xl = xl_^ + (6.20) 

where x{t) is the solution determined by the ODE (6.4). 

We want to hnd the solution to (6.2) of the form 

Sj = + Hi. (6.21) 

Then yi will satisfy the following equation: 


f Vi+i -yi = -Ei + e{G{-e^ Y!j=i + %•)) + 0{eH)) - G{-e‘^ Y!j=i E'( 7 (x°)))^, 

^ for i = 1, ■ ■ ■ , k — 1 

[ E •=! + 0(^2) E ■=! I%f = -E, + O(s^k), 

( 6 . 22 ) 

where 


Ei = -x°- eG(-e^ ^ H\7{x^))) 

i=i 


for i = 1, • • • , fc — 1, and 

k 

i=i 

First we show that the approximate solution we choose is indeed a good approximate 
solution, i.e. the error Ei is small enough. In fact, we have the following error estimate: 

Lemma 6.2. 

i 

Ei = - eG{-e^ ^ H'( 7 (x°))) = 0(e) (6.23) 

i=i 

for i = 1, ■ ■ ■ , k — 1, and 

k 

(6'24) 

i=i 

Moreover, the following estimate holds: 


k-l 

Ei®-I = c>(£). 

i=l 


(6.25) 



28 


W. AO, H. CHAN, AND J. WEI 


Proof. First for i = k — 1, we have 


fc-i 


xl - xl_^ - eG{-e'^ ^ i/'( 7 (x°))) 


t=i 


fc-i 


Pb) - eG{-e^ H'ilix^j))) 


i=i 


fc-i 




Since we choose the midpoint approximation, we have for i = 1, - ■ ■ ,k — 2, 


p{U+i) - = 0{h^), (6.26) 

i=i 


and 


k k—1 

YH'{l{x]))h = (YH'{l{x]))h-p{tk)) + {H'{'^{xl))h +p{tk)) 

i=i i=i 

= o(A^) + o(ihfM,fc = o(ihfM)A. 


In e 


Ine 


(6.27) 


By (6.26) and (6.27), and recall that pi, = 0{h), one can obtain that 


fc-i 

Ek-i = xl- xl_^ - eG{-e^ Y ^'( 7 ^°))) = 0{eh) 

i=i 


and 


Ek 


0(e2 


ln(— Ine) 
In £ 


Next by the equation satished by p{t), we can get that 


(6.28) 


p{ti) = 0(min{z, k — i + l}h) 


(6.29) 
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SO for z = 1, • • • ,k — 2 

rU+l 


X. 


J=1 


rn+i 1 p ^ 

I -n77c;(—p(())* - eG(-e^ 

Jti 


j=l 'I 

0(-^j^)(p(i.+i) - ^ H'(',{x]))h) + 


= 0 (— 


e/z 


-) + 0(e)(- 




min{z, A; — z + 1}^ min{z, A; — z + 1}' 


'min{z, A; — z + 1} 

= 0 (£). 

Moreover, from the above estimate, we have 

i 

= 0(e), for z = 1, • • • , A; — 1. 
i=i 

Finally, we will show that equation (6.22) is solvable. 


□ 


Lemma 6.3. There exists Sq > 0, such that for e < eo, there exists a solution {yi}i<i<k to 
(6.22) such that 

Iblloo < C'eln(-lne). 


(6.30) 


Proof. For \\y\\oo « ^1 In^l, we have 

i i 

eG{-e^ H'i-fix^j + yj)) + 0{eH)) - eG{-e‘^ ^ H'{-f{x°))) 
i=i i=i 


= _e(M!(2M)hL) + 0( 




-) + 0 (; 


e^z 


E;=i^'(7(:^°)) E;=i^'(7(:^?)) ^'(7(x°)) 

The equations (6.22) for z/j can be rewritten as follows: 

i) 


) 


_ E}^7^"(7(x°)fa ^ o(^Ms^) + 0( 


for z = 1, • • • , A: — 1 

E‘.i H"('i(-4))y, + Eh. = 0 (ek) + 




(6.31) 
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We will show that one can hrst solve 1/2 to t/k in terms of yi from the hrst k — 1 equations, 
and hnally solve yi by the k-th equation of (6.31). 

For 1 < I < io = {1 — 6) k where 5 > 0 is a small number to be determined later, we have 


yi+i -yi + ^Yl 


Ei=i H''{'j{x^j))yj 






ei 


i=l 




\y\'i<i + 




^e;=i^'(7(^°)) 


o(e) + V- 

' ' mi 


+ E^(— 

f ^ miT 




min{z, A; —-i + 1} ^min{z,/c — z + 1}' 

,el 


- 0(e) + 0(^)iy|i<i. 


where we denote by 
Moreover, 


\y\h<i<i 2 = sup \yi\. 


^^ T.UH"(^(x]))y, 

Y.]=iH'(-i(x^)) 


= £ 


5:o( 


Ay\ 


i<l 


i=l 


min{z, k — i + 1}' 


Thus one can get that for I < 

yi = yi + o(l)\y\i<io + o(l)\y\‘f<i^ + 0(e) 

So we can get that 

yi = (1 + o(l))yi + O(e)0 = 2, • • • ,io- 
For / > zq, we have the following: 

e;=i^^"(7(^°))%' 

yi+i -yi = -^ 2^ 


i=io+l 


e;=i^'(7(x°)) 


,eL 


+ 0(e) + + O(y)|?/|i<io + o(l)\y\i<ig 

= CoS\y\i^^i<i + 0 (| 2 /|iQ<i<i) + 0(\y\i<i2 + 0(e) 


(6.32) 

(6.33) 


for some Cq independent of £ and 6. So for io < i < k, 

yi = 0(yi) + C'o( 5 |z/|io<i<z + 0 (|z/|^g<j<;) + 0(e) 

If (5 > 0 is small such that CqS < then the above system is solvable with 

yi = 0 ( 1 / 1 ) + 0(£). 


(6.34) 


(6.35) 
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From the last equation, we have 

io k 

Y 1 H"h(xf))y, + 0(k\yA^) + 0(kE^) 

2=1 2 = 2 q -|-1 

= + o(l))?/i + 0{6k\yi\) + 0{ke) + 0{k\yi\‘^) 

2=1 

ln£ 

Thus by the assumption {Hi), the equation is reduced to 

2/1 = o{l)yi + 0{6)\yi\ + Odi/iH + C>(e ln(-Ine)) 

If we further choose S small enough but independent of e such that 0((5)|i/i| < \\yi\, it is 
easy to see that by contraction mapping, the above equation has a solution and satishes 


yi = 0(£ln(-lne)). 

Thus we get that there exists a solution to (6.22) with 

II 2 /II 00 < Ce ln(— In e) << e| Ine 
Thus we have proved the existence of solution to (6.22). 


(6.36) 


□ 


7. Appendix: Proof of Proposition 3.1 

In this appendix, we shall give a proof of Proposition 3.1. The proof is rather standard. 
It follows from arguments in [4] and [25]. It is based on Fredholm Alternative Theorem for 
compact operator and an a-priori estimates. 

First we need an estimate on the following matrix M dehned by 

Mij = / ZiZjdz, i,j = I,-■■ ,k. (7.1) 

Jne 

Lemma 7.1. For e sufficiently small, given any vector b G there exists a unique vector 
/5 G such that M/3 = b. Moreover, 

Halloo < C'llftiloo (7.2) 

for some constant C independent of e. 

Proof. To prove the existence, it is sufficient to prove the a priori estimate (7.2). Suppose 
that |/3j| = ||/3||oo, we have 

k 

2=1 
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For the entries Mj^, from the dehnition of A^, and the exponential decay property of Zj, we 
know that 

f f dvo 

Mu = / Z‘^dz = (1 + o(l)) / {^fdy > cq > 0, 

Jne 


and 

ElWI<C53e-—=0(1). 


Hence for e small, we have 

Coll^lloo < ColAI < ^ + \bi\ < o(l)||/3||oo + ll&lloo 

from which the desired resnlt follows. □ 


Second we need the following a priori estimate 

Lemma 7.2. Let h G with ||/i||* bounded and assume that ("0, {q}) is a solution to 

(3.5). Then there exist positive numbers Sq and C, such that for all e < Sq, for any points 
Pi, i = 1,... ,k given by (1.15) , one has 

mu<c\\h\u. ( 7 . 3 ) 

Proof. We argne by contradiction. Assnme there exist tf solntion to (3.5) and 

||h||* —0, IIV^II* = 1. 


We prove that 


Ci —0 for i = 1, • • • ,k. (7.4) 

Mnltiply the eqnation in (3.5) against Zj and integrate in we get 



L^ifZji^z) 



hZj + M{cj), 


By the exponentially decay of Zi, we hrst know that 


hZf < C\\hl, 
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Here and in what follows, C stands for a positive constant independent of e, as e ^ 0 . 
Secondly, by the equation satished by p.w{z — Pi), 


L^'ipZidz = 


[ (A?/> - i/; + p( V Ps 

Jn, 


heWiz-Pi))P ^'ip)Zidz 


» k 

/ (AZ, - Z, + p( V Pn, ^ w{z - Pi))^-^Zi)i;dz 

i=_k 


~ ^ - --pw{z-Pi)'P 


i=l 


dr 


dr 


< C 


'B\lne\{Pi) 


PJ I I O (e) + y; (2; _ p. )P 2 ^ y; (2; - P,. ) I I I 


d\w{z — Pi 


JOe/Rl ine|(Pi) -■ 

<CUUO{e) + 0{e^)) 

< CeWiplU 


J2 + 0(e) ^ e-^l"-OI] |V;|dz 


j¥^i 


i=i 


if we choose p small enough such that p — p > 2. This can be done since p > 2. 

Since M is invertible and ||M“^|| < C, we get that 

|Q|<C(||h||, + 0(e)||^||.). (7.5) 

Thus we get the validity of (7.4), since we are assuming |1'0|1* = 1 and ||/i||* 0. 

Let now /i G (0,1). It is easy to check that the function 

k 

W :=^ 

i=—k 


satishes 

in \ Uj=i_,,,^fcP(Pj, P) provided R is hxed large enough (independently of e). Hence the 
function W can be used as a barrier to prove the pointwise estimate 

I0l(^) < O ||i/’||L-(n(p,.,p)no,)j W{x ), (7.6) 

for all 2: G \ \JjB{Pj, R). 

Granted these preliminary estimates, the proof of the result goes by contradiction. Let us 
assume there exist a sequence of £ —)■ 0 and a sequence of solutions of (3.5) for which the 
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inequality is not true. The problem being linear, we can reduce to the case where we have a 
sequence tending to 0 and sequences such that 

and ||^/;W||* = 1. 


But (7.4) implies that we also have 



^ 0 . 


Then (7.6) implies that there exists such that 


(7.7) 

for some hxed constant (7 > 0. Using elliptic estimates together with Ascoli-Arzela’s theorem, 
we can hnd a sequence and we can extract, from the sequence — a subsequence 
which will converge (on compact) to i/joo a solution of 


{A-l+pw^ -000 = 0 , 


in M^, which is bounded by a constant times e ^^^1, with /i > 0. Moreover, since satishes 
the orthogonality conditions in (3.5), the limit function ipoo also satishes 



0 . 


But the solution w being non-degenerate, this implies that 0oo = 0, which is certainly in 
contradiction with (7.7) which implies that 0oo is not identically equal to 0. 


Having reached a contradiction, this completes the proof of the Lemma. 

□ 


We can now prove Proposition 3.1. 

Proof of Proposition 3.1. Consider the space 

U = {u e ■. [ uZi = 0, i = l,...,k}. 

Jfle 

Notice that the problem (3.5) in 0 gets re-written as 

0-|-A'(0) = h in H (7.8) 

where h is dehned by duality and iL : "H —)■ "H is a linear compact operator. Using Fredholm’s 
alternative, showing that equation (7.8) has a unique solution for each h is equivalent to 
showing that the equation has a unique solution for h = 0, which in turn follows from 
Proposition 7.2. The estimate (3.7) follows directly from Proposition 7.2. This concludes the 
proof of Proposition (3.1). 
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